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We report on a search for CP violation in the decay of and to if^Tr+Tr". The data come 
from an integrated luminosity of 9.0 fb""*^ of e^e~ collisions at -^5 « 10 GeV recorded with the 
CLEO II. V detector. The resonance substructure of this decay is well described by ten quasi-two- 
body decay channels plus a small non-resonant component. We observe no evidence for CP violation 
in the amplitudes that describe the decay K'g-K'^-n:~ . 

Phenomena that are not invariant with respect to charge conjugation and parity (CP) in strange and bottom 0, 
13 mesons are the motivation for numerous current and future experiments. Standard Model (SM) predictions for 
the rate of CP violation in charm mesons are as large as 0.1% for D'^ 7r+7r~7r° 0,0 but are considerably smaller, 

0(10"^), for D° — > KgTT^TT" where the dominant contribution is due to K° — K mixing 0. The Dalitz technique 0,0 
allows increased sensitivity to CP violation by probing the decay amplitude rather than the decay rate. Observation 
of CP violation in KgTT^Tr~ at current experimental sensitivity would be strong evidence for non-SM processes. 

The decay DK^ followed by a multibody decay, such as X^tt+tt", may elucidate the origin of CP 

violation in the B sector [lol |. 

We present the results of a search for CP violation in the amplitudes that contribute to 1?° — > Kg7r~^Tr~ . Previous 
searches for direct CP violation 0,0,0,0,0 in the neutral charm meson system set limits of a few percent. 
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This analysis uses an integrated luminosity of 9.0 fb~^ of e+e~ collisions at ^/s « lOGeV provided by the Cornell 
Electron-positron Storage Ring (CESR). The data were taken with the CLEO II. V detector T^ . 

The event selection is identical to that used in our previous analysis of K^tt^tt' - [13 which did not consider 
CP violation. We reconstruct candidates for the decay sequence I?*+^7r^D°, ^ if^7r+7r~. The charge of the 
slow pion (ttJ" or tt^) identifies the charm state as either or . Consideration of charge conjugation is implied 
throughout this paper, unless otherwise stated. 

We evaluate the energy released in the tt^D*' decay as (5 = M*— M— m^, where M* is the reconstructed mass 

of the ■k'^ K'g'K'^ t:~ system, M is the reconstructed mass of the K^t:'^t:~ system, and is the charged pion mass. 
We require the momentum to exceed 2.0 GeV/c. We reconstruct K^^tt^tt^ with the requirement that the 
daughter pion tracks form a common vertex, in three dimensions, with a confidence level > 10~^. Signal candidates 
pass the vertex requirement with 96% relative efficiency. Throughout this paper, relative efficiency is defined as the 
number of events in the data passing all requirements relative to the number of events when only the requirement 
under study is relaxed. 

Our silicon vertex detector provides precise measurement of charged tracks in three dimensions [T^ . We exploit 
the precision tracking by refitting the trajectory and tt^ tracks with a requirement that they form a common 
vertex in three dimensions. We use the trajectory of the K^'k^tt~ system and the position of the CESR luminous 
region to obtain the D° production point. We refit the tt^" track with a requirement that the trajectory intersect the 

production point. We require the confidence level of each refit exceed 10"''. The signal candidates pass the 
production and decay vertex requirement with 85% and 91% relative efficiency, respectively. 

We select 5,299 candidates within three standard deviations of the expected Q, M, and m^^o. We compute ctq, 
ctm and (Jm^^o from the trajectory reconstruction covariance matrices of the daughters of each £)*+ candidate. The 

distributions of Q and M for the and D° samples for our data are shown in Fig.^ We find 2,579 13° and 2,720 
candidates corresponding to an asymmetry of (— 2.7± 1.4 ± 0.8)%, a 1.7(7 effect. 
In Fig. 121 we plot rn^jzjg vs m^_|_^_ where m^5 denotes the 'right-sign' and corresponds to m^o^- for and mj^o^+ 

for Similarly, mws denotes the 'wrong-sign' and corresponds to m^o^+ for and m^o^- for D°. We study our 
efficiency with a GEANT (T^ based simulation of e+e cc events in our detector with a luminosity corresponding 
to more than three times our data sample. We observe that our selection introduces distortions due to inefficiencies 
near the edge of phase space, and fit the efficiency to a two dimensional cubic polynomial £{m\g, m^+^_). The 
reconstruction efficiencies for the and over the Dalitz plot are consistent with each other and we take them to 
be equal; so £{m\g, to^+^_) =£{m'^g, ml+^_) =£. 

Figure n shows that the background is small, but non- negligible, and we model our background as in Ref. hj. To 
model the background contribution in the Dalitz distribution we consider those events in the data and MC that are 
in sidebands five to ten standard deviations from the signal in Q and M and within three in mj^o. There are 235 

(579) 1?° and 210 (572) I?° candidates in this selection in the data (MC), about four times the amount of background 
we estimate from the signal region. We constrain the invariant mass of the sideband candidates to the I?" mass 
and compare with the background in our signal region from our simulation which also includes e+e~ annihilations 
producing the lighter quarks. We note that the background from the simulation is dominated by random combinations 
of unrelated tracks. The simulation predicts that the background uniformly populates the allowed phase space. We 
model this contribution to the Dalitz distribution by fitting the 13° mass constrained data sideband sample to a 
two dimensional cubic polynomial 'B(m^g,m^^^_ ). All parameters except the constant are consistent with zero as 
predicted by simulation, so B{mj^g,m'^^^^) = B(rnyyg,m'^^^_) = B. The normalization of the uniform background 
in the data exceeds the simulation by 21 ± 8%. Other possible contributions to the background, where tttt or Ktt 
resonances combined with random tracks fake a and real decays combined with random soft pions of the 
wrong charge. The latter, called mistags, are especially dangerous to our search for CP violation as candidates 
are misidentified as D^'s. Since only two of the Dalitz parameters are independent jl^l, mistags populate the Dalitz 
distribution in a known way that depends on the shape of signal - namely we interchange assignment of rriRs and 
rriws- When we analyze the Dalitz distribution, we allow a mistag fraction with an unconstrained contribution and 
we have looked for the contribution of a resonance, such as p° or iiT* (892) ~, plus random tracks to the background in 
the data and conclude that any such contributions are negligible. 

We parameterize the Kgir^n^ Dalitz distribution using the isobar model described in Ref. 0] where each resonance, 
j , has its own amplitude, aj , and phase, Sj . A second process, not necessarily of SM origin, could contribute to the j-th 
resonance. In general, we can express the amplitudes to the j-th quasi-two-body state as {aje^'^^^^'^^'>±bje^^^^^'^^'>)Aj — 
aje^''^''^'^'^l±^)Aj, with '+' for £>° and ' — ' for and Aj — Aj{m\g,rn^^) is the amplitude for resonance j as 
described in Ref. Thus aj and 5j are explicitly CP conserving amplitudes and phases, hj are explicity CP violating 
amplitudes normalized by the corresponding CP conserving amplitude Oj, and <j)j are explicitly CP violating phases. 
In the absence of CP violation bj and (f)j would be zero. The matrix elements M. and M. for the and D° samples. 




FIG. 1: Distribution of Q, (a) and (b), and M, (c) and (d), for the process Z?" and D° K%ti^-k . The candidates pass all 
selection criteria discussed in the text. 
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FIG. 2: Dalitz distribution for a) D° -» /sTg 7r"'"7r and b) D° — > Kgn'^n candidates passing all selection criteria discussed in 
Ref. [17]. The horizontal axis {mRs)^ corresponds to (mj^o^-)^ for D° and {'mjfO^+)^ for D°. 
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where ao and So parametrize the non-resonant amplitude, assumed to be CP conserving 
We perform an unbinned maximum likelihood fit which minimizes the function 
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define the normalization of the D", D", and background events. The signal fraction Fo and its error ap are determined 
from the fit to the combined and Z?" mass spectra, shown in Fig. QJ;) and d), to be 0.979 ± 0.015. The signal 
fraction F and the mis-tag fraction (1 — T) are consistent in the and samples and we take them to be equal. 
The signal fraction F, the mistag fraction (1 — T), and the parameters aj,Sj,bj, and that describe the matrix 
elements A4 and A4 are determined by the Dalitz plot fit. We test the performance of our fit by generating 100 
Monte Carlo samples from the results of our standard fit in Ref. The pull distributions of all fit parameters are 
consistent with unit Gaussian with zero mean indicating that the fit is not biased and that the errors are correctly 
computed. 

We begin our search for CP violation from the results of our standard fit in Ref. ^3 which clearly observed the ten 
modes, {K*-tt+, K*{U30)-tt+, i4'2*(1430)-7r+, X*(1680)-7r+, K°p, K°uj, ii:g/o(980), ii:g/2(1270), i^^/o(1370), and 
the "wrong sign" K*~^tt~) plus a small non-resonant component. First, we fit the and Z?" samples independently. 
The results of the and fits are consistent with each other and with our CP conserving result |T3|. Next, we fit 
the and samples simultaneously. This fit has 42 free parameters, ten CP conserving amplitudes and ten CP 
conserving phases which are the same in the and D'^ samples, ten CP violating amplitudes and ten CP violating 
phases which differ by a sign in the D'^ and D*^ samples, plus two normalizations for the combinatoric and mistag 
backgrounds. 

We report the CP conserving amplitude and phase, aj and Sj, in Table These results are consistent with our 
result in Ref. in which CP conservation was assumed. We report the fractional CP violating amplitude and CP 
violating phase, bj/aj and in Table ITU The three projections of the fit to the combined and samples and 
the difference between the D° and D° samples are shown in Fig. O We find the signal fraction F and mistag fraction 
(1 — T) to be, 97.1 ± 0.8% and O.O^Q g%, respectively. The confidence level of the fit, calculated directly from the 
likehhood function Hill, is 55%. 

The amplitude aj, on its own, is not a good estimator of the contribution of resonance aj to the total rate. The 
width of a resonance, interference with other resonances and the allowed phase must be considered. The fit fraction 
(FF), formulated to encapsulate all these effects, is commonly defined as the integral of a single amplitude squared 
over the Dalitz plot {mj^g vs rn^^) divided by the coherent sum of all amplitudes squared. We define the fit fraction 
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for the I?" and D'^ samples, respectively. The CP conserving and CP violating fit fractions are defined as the sum 
and difference of the numerators of Eq. |S1 and Eq. O respectively, divided by the sum of the denominators of Eq. |S1 
and Eq. 

The dominant constraint on CP violation is not due to limits on the CP violating amplitude squared but is due to 
the potential interference of a CP violating amplitude with a well determined CP conserving amplitude. We define 
the CP violating interference fraction (IF) as 

_ l/Efc i^ake'^" sin(0fc + 0j)A) bjAjdm\sdml^\ ^^^^ 
(/ \M\^ dm\sdml^ + / l^f dm\sdml^^ 

where for the non-resonant component = 1. The value of bj determined by our fit is constrained by terms in the 
likelihood function proportional to \bj\^ and ake^^^''bj which are sensitive to both CP violation in the direct decay 
to a given submode and possible CP violation in interference with other modes, respectively. The CP violating fit 
fraction defined by Eqs. |H1 and |51 is sensitive to CP violation in decay. The CP violating interference fractions of 
Eg. IIUI sum over the contribution proportional to ake^^^''bj so are sensitive to CP violation in interference. The 
phases are important and allow the possibility of cancellation in this sum. This makes the IF a better representation 
of the impact of CP violation on the rate of decay than the CP violating FF. To quantify a fractional CP violation 
each decay channel, we define 

AcP, =IF,/FF,, (11) 



where FFj are the CP conserving fit fractions. 
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Note that the observation of CP violation with a significantly non-zero bj/aj or does not imply that any of the 
derived quantities, the CP violating FFj, IFj or Acp^ will be non-zero. In particular, the CP violating effects that 
are locally non-zero can integrate to zero, thus IFj and Acp^ are more sensitive measures of CP violation than the 
CP violating fit fraction. 

We use the full covariance matrix i22j from the fits to determine the errors on fit fractions and the CP violating 
interference fractions so that the assigned errors will properly include the correlated components of the errors on the 
amplitudes and phases. After each fit, the covariance matrix and final parameter values are used to generate 500 
sample parameter sets. For each set, the fit fractions are calculated and recorded in histograms. The statistical error 
on the fit fractions is then extracted from the histograms. In Tabled we report the results for the CP conserving fit 
fractions, and the 95% upper limit for CP violating contributions are given in Table UTI The fit fractions for the D'~' 
and samples are given in Table Hill An alternative measure of the rate of CP violation in a given submode is the 
asymmetry between the and fit fractions, which arc also given in Table III. The "Fit Fraction Asymmetry" is 
similar in sensitivity to Aqp defined by Eq. If 11 

The common evaluation of the integrated CP asymmetry between normalized amplitudes squared across the Dalitz 
plot is sensitive to an asymmetry in shape between the and 13° samples and is defined as 



We obtain Acp = — 0.009±0.02f _QQ4!jrQQ37 where the errors are statistical, experimental systematic and modeling 
systematic, respectively. 

We consider systematic uncertainties from experimental sources and from the decay model separately. Our general 
procedure is to change some aspect of our fit and interpret the change in the values of the amplitudes, phases and 
fit fractions in the non-standard fit relative to our nominal fit as an estimate of the systematic uncertainty. The 
impact of systematic uncertainties on the upper limit of the interference fraction, CP violating fit fraction and Acp 
reported in Table ITU are estimated by recomputing the statistical error on these quantities with the covariance matrix 
of the non-standard fits using procedure described above. Contributions to the experimental systematic uncertainties 
arise from our model of the background, the efficiency, and biases due to experimental resolution. The background 
is modeled with a two dimensional cubic polynomial and the covariance matrix of the polynomial coefficents, both 
determined from a sideband. Our nominal fit fixes the coefficients of the background polynomial, and to estimate the 
systematic uncertainty on this background shape we perform a fit with the coefficients allowed to fioat constrained by 
the covariance matrix. Similarly we perform a fit with a uniform efficiency, rather than the nearly uniform efficiency 
determined from the simulation, as estimates of the systematic uncertainty due to the efficiency. We also perform 
fits where the background normalization and the efficiency of the and D*^ samples are determined separately. We 
compute the overall normalization by evaluating the integrals in Eqs.ElandElusing Gaussian quadrature to interpolate 
between points on a finite grid across the Dalitz plot. To study the effect of the finite resolution our experiment has 
on the variables in the Dalitz plots we vary the granularity of the grid used to compute the overall normalization. 

We change selection criteria in the analysis to test whether our simulation properly models the efficiency. We 
introduce a track momentum cut of 350 MeV/c to avoid the difficulty of modeling our low momentum tracking 
efficiency. We expand the signal region from three to six standard deviations in Q, M, and TO/^o, and increase the 
Pd- cut from 2.0 to 3.0 GeV/c. These variations to the nominal fit are the largest contribution to our experimental 
systematic errors. 

Contributions to the theoretical systematic uncertainties arise from our choices for the decay model for —> 
KgTr~^TT~ . The standard value for the radius parameter for the intermediate resonances and for the is 0.3 fm 
and 1 fm, respectively. We vary the radius parameter between zero and twice the standard value. Additionally, we 
allow the masses and widths for the intermediate resonances to vary within their known errors pol l24l l25j | . 

We consider the uncertainty arising from our choice of resonances included in the fit. We compared the result of 
our nominal fit to a series of fits where each of the resonances, a or /o(600), /o(f500), /o(f7fO), p(I450), p(f700) 
were included one at a time. We also considered a fit including both the /o(600) and /o(I500) resonances. These 
variations to the nominal fit result in highly asymmetric variation in fit parameters and are the largest contribution 
to our modeling systematic error. 

We take the maximum variation of the amplitudes, phases and fit fractions from the nominal result compared to 
the results in this series of fits as a measure of the experimental systematic and modeling systematic uncertainty. 

In conclusion, we have analyzed the resonant substructure of the decay or KgTT^TT~ using the Dalitz-plot 

analysis technique and searched for CP violation in the amplitudes and phases of the ten clearly observed intermediate 
resonances. Our results, shown in Table UTI and Table ITTTI are consistent with the absence of CP violation. We find 
the CP asymmetry in the fit fractions for each decay channel to be in the range < (3.5 to 28.4) x 10^^ at the 95% 
confidence level. We find the CP asymmetry in the interference fractions for each decay channel to be in the range 




(12) 
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FIG. 3: Projections of the results of the fit described in the text to the Kg'K'^n~ Dahtz distribution showing both the fit 
(histogram) and the data (points). The combined -D° and D° samples are shown in (a), (b) and (c). The difference between 
the D° and D° samples is shown in (d), (e) and (f). 
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< (0.4 to 22) X 10"'^ at the 95% confidence level. We find the ratio of the CP violating to CP conserving rate for each 
decay channel, to be in the range < (0.3 to 150)% at the 95% confidence level. We find Acp which is the asymmetry 
between normalized squared amplitudes integrated over the entire Dalitz plot to be -0.009±0.021+°;°4!^°;°3^. 
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TABLE I: CP Conserving Parameters. Errors are statistical, experimental systematic and modeling systematic, respectively. 
The CP conserving fit fraction is computed from Eqs. |5] and |3] following the prescription described in the text. 
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U.04ztU.UZ Q Q2 UlU.Ul 


iO5=t0_4_4 


A 9-1-0 f^^-^-l-OO 


Z°/2(1270) xB(/2(1270)^7^^ 




o.79±o.2SO:?jlg:^,? 


qnfi+13+7 +54 
■J>-'U_j^g_23_3 


0.36±0.22t°:?^«:?« 


Z°/o(1370) X B(/o(1370) 7r^ 




1 74+n 1 ■^0.23+0.10 
1. Cy:=tU.li_Q 22-0.20 


85±5t5ti2 


9.8±l.#^2li2 9 


ii:o*(1430)-7r+ X 5(7^0* (1430)" 


^ k\-) 


1 qq-l-n 1 1+0 04+0.30 


■|_LC+&f8 
i-a:o_4_l4 


7.2±o.^;:?tj:^ 


Js:2*(1430)"7r+ X B(/C2*(1430)" 




n Q/1+0. 12+0. 13+0.00 
^■''^-0.11-0.10-0.06 


qqir+9+0+4 
•JOO_g_g_23 


1 1+0 2<-0-3+o.4 


Ji'*(1680)-7r+ X B(A'*(1680)- 




r. 4qf0.67+0.74 , o q 
O-*y_0.65-0.94='='J-y 


175±^iilf 




K tt'^tt~ non-resonant 




r, Qo+0. 34+0. 59+0. 92 
'5-0. 31-0. 17-0. 53 


q/iq+23+58+74 

'J^'J- 16-23-17 





TABLE II: CP Violating Parameters. Errors are statistical, experimental systematic and modeling systematic, respectively. 
The interference fraction, CP violating fit fraction and Acp, computed from Eg.lStEg. lllI following the prescription described 
in the text, include statistical and systematic efi'ects. 



Component 



A'*(892)+7r- X B(A'*(892)+ ^ K''-k+) 

K^uj X B{ijj tt+tt") 

if*(892)-7T+ X B{K*{892)- -^Tfn-) 

K°fo{980) X B(/o(980) ^ tt+tt") 

ZV2(1270) X B(/2(1270) ^ TT + TT") 

ZVo(1370) X B(/o(1370) ^ tt+tt") 
A'o(1430)-TT+ X B (/Co* (1430) - 
ii'J(1430)-TT+ X S(74'2*(1430)~ 
A'*(1680)"TT+ X B(A'*(1680)" 



Ratio {bj/aj 



Phase (03 ) (° 



— 1 9fu.2i+u.oy+u.ii 

'-'•1^-0.22-0.15-0.03 

o.ooi±o.o2i°:S9t;5:E]?? 

_n 1 /pi-o. 10+0. 11+0. 01 

'-'•1^-0.11-0.01-0.02 

^Om2±OMtlZt'-°oZl 

-o.o4±o.ofi°:jn°° 

Q j^gf 0.28+0.15+0.11 



Interference 
Fraction 

(95% 



CP Violating 
Fit Fraction 
Upper Limits) 



K TT 
K TT 
K TT 



o.os!;o 



0.27-0.37-0.18 
0.06+0.01+0.06 
05-0.11-0.03 



-0 02+0 06^0.04+0.00 

U.UZZ1ZU.UD_Q Q2_o.oi 
-0 O^^l + O 12^0.04+0.04 
0.28+0.05+0.02 



-O.2OI; 



0.27-0.22-0.01 



-+21+13+18 
"-22-35-4 


< 


3.0 


X 


10" 


-3 


< 7.8 X 10~ 


4 


< 


92 


1+16+9 +21 
1-18-31-3 


< 


0.7 


X 


10" 


-3 


< 4.8 X 10" 


4 


< 


0.3 


0+17+8 +20 
°-19-30-3 


< 


0.4 


X 


10" 


-3 


< 9.2 X 10" 


4 


< 


4.5 


q+16+9 +21 

•^-18-30-3 


< 


2.1 


X 


10" 


-3 


< 3.5 X 10" 


4 


< 


0.3 


gf 16+10+20 
^-17-29-3 


< 


4.2 


X 


10" 


-3 


< 6.8 X 10" 


4 




10.4 


9r+19+12+20 
^^-20-32-2 


< 


4.4 


X 


10" 


-3 


< 13.5 X 10 


-4 


< 


150 


0+15+10+20 
"-17-28-4 


< 


22 


X 


10^ 


3 


< 25.5 X 10 


~4 


< 


21 


q+17+13+23 
■^-19-36-2 


< 


9.1 


X 


10" 


-3 


< 9.0 X 10" 


4 


< 


14 


q+17+10+21 
■^-18-31-2 


< 


2.2 


X 


10" 


-3 


< 6.5 X 10" 


4 


< 


22 


q+19+20+27 
•J- 20- 25- 2 


< 


19 


X 


10" 


3 


< 28.4 X 10 


-4 


< 


92 



TABLE III: D° and Samples. Errors are statistical, experimental systematic and modeling systematic, respectively. The 
fit fraction is computed from Eqs. |H| and 1^1 following the prescription described in the text. The Fit Fraction Asymmetry is 
computed as the difference between the D° and Fit Fractions divided by the sum. 



Component 


D° Fit Fraction (%) 


D° Fit Fraction (%) 


Fit Fraction Asymmetry (%) 


A'*(892)+TT" X B(A'*(892)+ ~* 




-) 


27±0 20*'"'^^+"-^" 
u.zrzizu.zu_o ii_o.o3 


41 + 21+"-^'^+"-"y 
u.y:i=cu.zi_o 00-0.05 


—21 +4.2^'i''+2^ 

^l=cy:^_28-4 








27.5±l.ftJ-^l°-^7 


25.9+l.^°1+2 7 


3.i±3.^?-i«:^ 


K^iu X B{iu tt+tt") 






61±0 2#0-32+o.i5 

U.Oi=l=U.zq:_o 09-0.09 


i.o3+o.3i«:rg:J? 


-26+2#fl4 


ii'*(892)-TT+ X B{K*{892)- 




-) 


68.0±1.8+^-^+^-^ 


64.7±1.^°J«;« 


2 5-1-1 qfi.5+2.9 

z.on: j..t;_o y_o.3 


K°fo{980) X B(/o(980) -> tt+tt 






4.0±0.^J-+°-i 


4.4+0.7li:3±0.2 


-4.7+ll.Ct?*4^1°j 


;^°/2(1270) X B(/2(1270) ^ TT" 






0.49+0.411^2^12:^^ 


24+0 2'^<-0-''i+0 06 
u.z:q;zizu.zo_o 03-0.09 


34+5ll7^l34 


ZVo(1370) X B(/o(1370) tt" 






ii.7+i.gtl^2l°2:3 


8.2±i.tlitl:l 


18±l(tl{tf 


7^0 ( 1430) "TT+ X 5(7^0* (1430)" 






7 111 1+0.9+1.5 
r.lzlzi.i_o 7_o.6 


7 9-1-1 1+0.4+1.2 
' •^=^^1-1-1.4-0.7 


-o.2±n.3^li+\-l 


ii'|(1430)"TT+ X B(A'2*(1430)" 




^") 


1.0+0.4+0.#°-5 


1 1+0 3+0.5+0.4 


-7±2tl,t\' 


A'*(1680)-TT+ X B(7C*(1680)" 








3.2±0.8^°-J+1i 


-36±19^l+l 



